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The moving-mirror problem is microscopically formulated without invoking the external boundary 
conditions. The moving mirrors are described by the quantized matter field interacting with the 
photon field, forming dynamical cavity polaritons: photons in the cavity are dressed by electrons in 
the moving mirrors. The effective Hamiltonian for the polariton is derived, and corrections to the 
results based on the external boundary conditions are discussed. 
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I. INTRODUCTION 

When mirrors are closely placed, the attractive or 
repulsive force between them is observed. This phe- 
nomenon, known as the Casimir effect is explained 
by the fact that the vacuum state of the electromagnetic 
(EM) field in the presence of the mirrors is modified from 
that of the free space, and the vacuum fluctuation energy 
depends on the positions of the mirrors. On the other 
hand, when the mirrors move very rapidly, quantum state 
of the EM field cannot adiabatically follow the instanta- 
neous vacuum state for each position of the mirrors, re- 
sulting in the creation of photons. Such excitation of the 
quantum field caused by non-adiabatic change of the vac- 
uum state B |[ ^] is referred to as the dynamical Casimir 
effect (DCE), and there have been numerous investiga- 
tions into this subjec t [f| 



tions into tms subject \f 




e.g., spectral properties of created photonfTH, radiation 
pressure on a moving mirror ||, [?], ^| , squeezing in the ra- 
diation field H |To), effective Hamiltonian approach lO, 



|12| , |13| , |14j, time- varying refractive index [15, [16 
radiation from moving dielectrics | 19|, 20|, 21 



17 

23 



influence of finite temperature p5. |26rp7|, and relation 
with sonoluminescence |28|, p9[ . 

In the most of the previous works, the moving mir- 
rors have been treated as the moving boundary condi- 
tions such that the transverse components of the electric 
field operator vanish at the mirror surfaces in their rest 
frames. Such external boundary conditions, of course, vi- 
olate the commutation relation of the EM field operators 
at the boundaries, and the incarceration of the photon 
field between the moving mirrors causes the temporal 
change of the Hilbert space. Thus, the 'classical' exter- 
nal boundary conditions involve quantum mechanical im- 
perfections. To circumvent these conceptual difficulties, 
the various results of the moving-mirror problem based 
on the boundary conditions should be examined and de- 
rived as some limiting case of more elaborate models. 
Several studies have been done towards this direction 
by considering moving matter with finite refractive in- 
23 , p3[: formulation of the problem 



dex |9j, J2C|, ^ , £2] u 

and radiation spectrum [p_9| , radiative reaction on the 



in dielectrics p4| . 

Recently, Koashi and Ueda ^] formulated the static 
Casimir effect based on a combined system of the EM and 
matter fields, and showed that the both fields participate 
in the vacuum fluctuations inducing the Casimir force. 
Although the quantum theory of the systems in which 
the EM field and matter interact with each other has 
been developed by many authors 32, [33. 34, pG| p7[, 
such an approach to the DCE has hardly been made so 
far 0. 

The aim of the present paper is to formulate the DCE 
in a moving-matter system in terms of the quantized 
field-matter theory. The EM field attenuates by coupling 
with the matter field inside the mirrors, and therefore 
no external boundary conditions are required. In other 
words, the EM field confined in the resonator is dressed 
by the matter field inside the mirrors, forming cavity po- 
laritons. Non-adiabatic movement of the mirrors excites 
the cavity polaritons and this phenomenon may be called 
the DCE of polaritons. In this paper, we derive the effec- 
tive Hamiltonian for polaritons, apply it to the moving- 
mirror problem, and compare the result with that based 
on the external boundary conditions |]l2| . 

This paper is organized as follows. In Sec. II, the 
moving- mirror problem is briefly reviewed. In Sec. Ill, 
we formulate quantum theory of field-matter interact- 
ing systems, in which matter is allowed to move. In 
Sec. IV, we derive the effective Hamiltonian for polari- 



19, E> 



in two and three dimension 



the dispersive mirror |20 



21L radiation 



density variation 



tons in the moving-mirror system, and apply it to the 
one-dimensional case. Final section presents the sum- 
mary of this paper, and some complicated algebraic ma- 
nipulations are relegated to appendices. 



II. BRIEF REVIEW OF THE 
MOVING-MIRROR PROBLEM 

We briefly review the moving-mirror problem || to 
make the present paper self-contained and to fix the no- 
tation. The simplest system consists of two perfectly 
reflecting mirror plates placed in parallel as illustrated in 
Fig. |l|. The mirror at the origin z — is fixed and the 
other at the position z = L(t) is allowed to move. The 
system is assumed to be uniform in the x and y direc- 
tions, and we consider only one component of the vector 
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FIG. 1: Schematic illustration of the one-dimensional 
moving-mirror problem. The left mirror is fixed at z = 0, 
and the right mirror moves along the z axis. 



potential, say the x component A x (z,t), without loss of 
generality. The vector potential in the Coulomb gauge 
obeys the wave equation as (we omit the subscript x of 
A x from now on) 



d 2 A(z,t) d 2 A(z,t) 



2 at 2 



dz 2 



(1) 



and the boundary conditions are imposed as A(0, t) = 
A(L(t),t) = 0, which guarantee that the transverse com- 
ponents of the electric field vanish at the surfaces of the 
mirrors in their rest frames. The field operator of the 
vector potential A(z,t) in the Hcisenberg representation 
can be expanded as 



A(z,t) = W— [a n {t)f n (z,t) + a+ (*)/„* (z,t)] , (2) 
„ V 2£ ° 

where d n and ajj are the annihilation and creation oper- 
ators of photons of the nth mode. 

One approach Q to this problem is to fix a n and aj, 
to ones at t = 0, and evolve the function f n (z, t) as 



d 2 f n (z,t) d 2 f n (z,t) 



c 2 dt 2 



dz 2 



(3) 



with the boundary conditions 

fn(0,t)=fn(L(t),t) = 0, 



(4) 



which ensure that A(z, t) obeys the wave equation (Q) 
and the boundary conditions A(0,t) = A(L(t),t) = 0. 
When L is constant, the function f n {z, t) is given by 



fn(z,t) 



aj n L 



' sin k n z, 



(5) 



where k n = nn / L and oj n = ck n . When the characteristic 
time of the mirror motion is much larger than L(t)/c, 
f n (z,t) adiabatically follows the mode function for each 
L(t) as 



fn(z,t) 



u n (t)L{t) 



3- J J>" (r)dT sin£:„(i)z, (6) 



where k n (t) = nir/L(t) and uj n (t) = ck n (t). When the 
mirror moves much faster, the adiabatic theorem breaks 
down, and f n [z, t) evolves in a more complicated manner. 
The time evolution of the system in this approach is thus 
not generated by a predetermined Hamiltonian but by 
the classical equation of motion (^) and the boundary 
conditions (|j), by which the time evolution of A(z, t) in 
the Heisenberg representation is obtained. 

Another approach to the moving-mirror problem is the 
method of an effective Hamiltonian 1^, [l4|, where 
the function f n (z,t) is fixed to the mode function for 
each L(t) as 



fn(z,t) 



w„(t)L(t) 



sin k n {t)z, 



(7) 



and a n {t) and ajj(i) are time dependent. This time evo- 
lution is described by the Hcisenberg equation as 



ih^ = [a(t),Hf(t)], 



(8) 



where the effective Hamiltonian for an operator O is de- 
noted by Hg(t). In the one-dimensional moving-mirror 
problem, the effective Hamiltonian for a n and ajj is ob- 
tained as [12] (see Appendix |X] for the derivation) 



If the first term is dominant on the right-hand 
side of Eq. (^), i.e., L/L <C w n for any n, 
the time evolution operator approximately becomes 



exp 
orem 



, showing the adiabatic the- 



One of the advantages of the effective Hamiltonian ap- 



3 



proach is that we can study the state evolution in the 
Schrodinger representation by 



^§- t \m) 



Hf(t)\m). 



(10) 



In this case, we note that the operators a n and fit can be 
interpreted as the annihilation and creation operators of 
photons of the nth mode defined in the interval < z < 
L(t). Another advantage of the effective Hamiltonian 
approach is that we can understand easily what kinds 
of elementary processes occur. For example, the second 
term of the effective Hamiltonian ([)]) has the form of the 
parametric process, suggesting that the squeezed state is 
produced when the mirror oscillates at frequency 2uj n . 
The third term induces pair creation and annihilation of 
photons and energy transfer between different modes. 

In the above approaches, the boundary conditions that 
the filed operator A(z,t) vanishes at z = and z = 
L(t) are imposed. Because of these boundary conditions, 
the canonical commutation relation [A(r, t), E(r', t)] — 
— iK/eoSrir— r') does not hold at the boundaries, where 
8t is the transverse delta function. The field operator (j^) 
and the Hilbert space on which it operates are defined 
only within the interval < z < L(t), and then the 
Hilbert space varies accordingly. 



III. FORMULATION OF THE FIELD-MATTER 
INTERACTING SYSTEMS 

A. Field representation of systems 

We start from a classical microscopic model, in which 
polarizable atoms in the matter interact with the EM 
field. We suppose that the ith atom consists of an elec- 
tron with charge — e and mass m e at the position and 
an ion with charge +e and mass M at the position R,. 
The relative vector between the electron and the ion is 
denoted by = — Ri, and the center-of-mass vector 
is Hi = (MRj + m e Yi)/{M + m e ) ~ Ri. We assume 
that the center-of-mass vectors S-i(t) are given functions 
of time when the matter is moved. The kinetic energy 
of the zth atoms is expressed as m e r 2 /2 + MRf/2 = 



(M+m e )Sf/2+mxf/2, where (M- 



/2 is a known 



function of time. Therefore, we consider mxf /2 alone in 
the dynamics, where m = Mm e /{M + m e ) ~ m e is the 
reduced mass. Furthermore the electrons and ions are as- 
sumed to be bounded by the effective potential mf2 2 x 2 /2. 
The Lagrangian for this system is then given by 



L 



dr 

E 



|E 2 (r,t)--L B 2 (r,t) 

m . 2 mfl 2 2 
~2 Xi ~ 



■ J2 [Ri- A(Ri,t)-Ti- A(ri,t)] 



<f>{n,t)] 
(ii) 



where E = -V</> - dA/dt and B = V x A. The Euler- 
Lagrange equations are obtained as 



e V ■ E(r, t)~eJ2 - Ri) - 5(r - r,)] 

i 

dE(r,t) 1 „, x 

£o p.' ' VxB r,t 

at fi 

+e ^ [BhS(t - Ri) - r*<5(r - r<)] = 0, 

i 

Xi = -0 2 Xi - — [E(r 4 , t)+r l x B(r,, t)} 



e 

M 



E(R{,t) + R t x B(Ri,t) 



0, 

(12a) 
(12b) 

(12c) 



The first and second equations are the Maxwell equa- 
tions, and the third one describes the motion of charged 
particles in the EM field. 

We rewrite the above particle picture of the polar- 
izable atoms in terms of the field picture. When the 
difference between adjacent polarizations |xi+i — Xj| is 
much smaller than the characteristic amplitudes of |x«| 
and |xj+i|, namely the polarizations change smoothly in 
the lattice scale, we can replace Xi with the polarization 
field X as 

Xj(t)=>.X(r,t). (13) 
The density of the polarizable atoms is replaced as 



^5(r-H,(t))=^p(r,i). 



(14) 



The polarization field X and the density p vanish outside 
the matter. The velocity of the matter &i(t) is denoted 
by v(r,t), which is defined only inside the matter. The 
time derivative Xj(i) should be replaced by dX.(r,t)/dt 
with 



d 
~dt 



d_ 
dt 



v(r,t) ■ V. 



(15) 



In the present paper, for simplicity, we consider the case 
in which the matter is allowed to undergo only transla- 
tional motion. Rotations and deformations of the matter 
complicate extremely the formulation and are not consid- 
ered. In this case, the velocity of the matter is uniform 
in each object, i.e., Vv(r, t) = 0. For example, in the 
case of Fig. |l|, v(z,t) = for z < and v(z, t) = L(t) 
for z > L(t). The density p(r,t) becomes a function of 
r — vt with v the velocity of each object, giving 

*M-(| +T .v),<„i,-o. m 

Using Eqs . (H^ ) and (|T4|), the first summation in the 
Lagrangian (|ll|) is replaced by 



drp(r, t) 



<£X(r, t) 
dt 



n 2 X. 2 (r,t) 



(17) 
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The interaction terms in the Lagrangian (|llj) can be 
rewritten as 



/ m \ 
r Xj, i 



r H Xj,t 



R A (r- ^jj x ^) + ' A ( r - 



-Xi,t 



xS(r- Si) 



-x, • V0(r,t) +x, ; • A(r,t) 



-Si • [ Xi - V]A(r,t) 



*(r-S,), 



(18) 



where in the last line we assumed that cj> and A are slowly 
varying functions in the scale of |xj|, and ignored the sec- 
ond and higher order of |xj|. Applying the replacements 
© and (|l|) to Eq. © yields 

J drp(r,<)|-X(r,t) • V0(r,t) + • A(r,t) 

+v(r,t) ■ [X(r,t) • V]A(r,t) 
= J drp(r, t)X(r, t) ■ [E(r, t) + v(r, t) x B(r, t)] + 

j t J drp(r,t)X(r,t) ■ A(r,t), (19) 

where the second integral in the second line can be ig- 
nored because the total derivative term in the La gran gian 
is irrelevant in the dynamics. Thus, from Eqs. (|17|) and 
@, the La grangian for the system is obtained as 



L 



m , 



( dX(r,t) 
{ dt 



fl 2 X 2 (r,t) 



-ep(r,t)X(r,<) • [E(r, t) + v(r, t) x B(r,t)]|20) 

We note that this Lagrangian reduces to the one used in 
the static and uniform dielectrics |3(|, when the matter 
is fixed and the density p is uniform. 

The conjugate momenta of A and X are given by 

5L 

n(l%t) = 8[d t A(r,t)] 

= -eoE(r,t)+ep(r,t)X(r,t), (21) 
,r/ s ^ , ^X(r,f) 

The Euler-Lagrange equation for 4> reads 

V-D(r,*)=0, (23) 



where D = —II can be regarded as the electric displace- 
ment with polarization — epX. Adopting the Coulomb 
gauge (V ■ A = 0) in Eq. fl23]), we can write the electric 
potential </> as 

0(r,f) = -f=^V-[p(r,t)X(r,t)], (24) 

where ^/(r) = -/ <fr7(r')/(47r|r - r'|). The Euler- 
Lagrange equations for A and X are obtained as 



£o- 



fl E(r,t) 
dt 



d 2 X(r,i) 
1 eft 2 



1 s , ,dX(r,i) 

= — V x B(r,t) + ep(r,t) ^' ; 

-ev(r,i)V.|>(r,*)X(r,i)], (25) 

= -mft 2 X(r,<) 

-e[E(r,t)+v(r,t)xB(r,| (26) 



Equation ([25|) corresponds to the Maxwell equation V x 
B//io = J + dD/dt with current J = evV ■ (pX) — e(v • 
V)pX, which satisfies V • J = 0. Equation (Eq) describes 
the dynamics of polarization moving in the EM field. Us- 
ing Eqs. (E0h-(E3|), we obtain the Hamiltonian as 



L 



H = dr 



(9/; 



/*{^P(r,t) 



1 



ep(r,t)X(r,t)]- + — B 2 (r,t) 



:Y 2 (r,t) 



2mp(r,t) v ' 2 
+ep(r, t)X(r, t) ■ [v(r,f) x B(r,i)] 

-Y(r,t) • [v(r,i) ■ V]X(r,i) L. 



p(r,t)X 2 (r,t) 



(27) 



Equations (|2l|), ([22|), ([25|), and ([26|) are derived as the 
canonical equations of this Hamiltonian. 

We note that the Hamiltonian ( p7| ) reduces to the 
macroscopic model described in terms of the dielectric 
constant e Q §l], |§, |1| ||, |j| , when the dynamics 
of polarization can be eliminated. The Hamiltonian in 
this model is given by 



H = J dr± (E • D 



4 



B H) 



1 



E D H B 

Mo 



1 - ~ ] v i D B i 



where we used the relations for moving medium J3 

D = e~E+ (e-e )v x B + 0(v 2 ), 
B = p H + (e-eo)E x v + 0(t> 2 ). 

If we neglect the left-hand side of Eq. p6|) , we get 



mQ 2 



(E 



x B) 



(28) 



(29a) 
(29b) 



(30) 
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This relation and D = e E — epX yield Eq. (29a) with 
e = l+e 2 p/s mSl. 2 . Substituting Eq. (30) into our Hamil- 
tonian (27), and dropping the terms including Y, we 
obtain Eq. (|28|). Our model, therefore, reduces to the 
above dielectric model, when the dynamics of polariza- 
tion is neglected. This corresponds to the case in which 
the relevant frequencies of the EM field are much smaller 
than the characteristic frequencies of matter Q and lu p . 



B. Quantization: polaritons 

We find out the normal mode of the field-matter cou- 
pled equations, in which the positions of the matter are 
fixed, i.e., v = 0. We denote this fixed matter con- 
figuration as M, and the mode functions and frequen- 
cies depend on it: A n (r,M), ui n (Ai), where n is 
the index of the mode. For brevity, we omit the argu- 
ment M below. Substitution of (h/2ij n ) 1 / 2 A n (r)e~ iWnt , 
i(?i W „/2) 1 /2n„(r) e —" t , z(^„/2) 1 /2 X „(r)e-^ t , and 
(n/2w n ) 1/2 Y„(r)e-^" t into A, IT, X, and Y in Eqs. ©, 
(p2|), (p5|), and ( |2"6| ) for v = yields the normal- mode 
equations 



II n (r) = -e A„(r) 



+< (l -V^V-) [/'!r)X„ (!•)]. (.Tlal 



^n n (r) = 

Y„(r) = 
Y„(r) = 



— V 2 A„(r), 
Mo 

mw 2 p(r)X„(r), 

e 2 

mfl 2 p(r) H p 2 (r) 

£ o 

--p(r)n„(r), 



X„(r) 



(31b) 
(31c) 

(31d) 



where we used Eq. (|24|). We take A„, n n , X„, and Y„ to 
be real without loss of generality. These mode functions 
can be shown to satisfy the orthonormal relation (see 
Appendix |Bj) 

J dr [A„(r) • II„,(r) - Y„(r) • X„,(r)] = -<W- (32) 

In terms of the mode functions, we can expand the 
time evolution of the fields as 



A(r, t) =J2\I ^ A "^ ( & « e """ t + c - c . ( 33s 



n(r, t) =J2 W^r n »( r ) ( 6 « e """* - c - c -) >( 33b ) 



X(r,t) = ^iJ^X n (r)(6 n e- iw « t -c.c.),(33c) 



Y(r,i) = ]T \/^7 Y »( r ) (p n e- iUnt + cc.) , (33d) 



where b n is the complex amplitude of each mode, and cc. 
indicates the complex conjugate of the previous term. 
Substituting the expansions (|3^) into the Hamiltonian 
(P7|) for v = 0, denoted by H v=0 , and using Eqs. (|3l]), 
we obtain 



H v =o 



[ hu} n b*b n . 



(34) 



Following the standard quantization procedure, we re- 
place the c-numbers b n and 6* with the Bose operators 
b n and 6^ satisfying the commutation relation [b n , b n ,\ = 
S nn ' ■ The elementary excitations created by w n can be re- 
garded as polaritons, since they are linear combinations 
of the photon and polarization fields. We should note 
that the operators b n also depend on the matter configu- 
ration Af, and are to be expressed as b n (A4) in full detail. 
The field operators in the Schrodinger representation can 
be expanded as 



A w = E\/^- A «( r )( s «+^)> ( 35a ) 



n(r) = £n/-^II n (r)(S n -St), (35b) 



X(r) = ^W^X„( r )(&„-&t), (35c) 



Y w = Ev^r Y «( r )(^+^)- ( 35d ) 



In the above argument, we assumed the discrete spec- 
trum of polaritons. When the spectrum is continuous, 
the continuous index of the mode, such as wave number 
k, should be used instead of n, and the summation J2 n 
should be replaced by an appropriate integral. 



IV. FIELD-MATTER FORMALISM OF THE 
MOVING-MIRROR PROBLEM 

A. The effective Hamiltonian 

In the previous section, the polaritons were derived for 
fixed matter configuration M . The number states of the 
polaritons are eigenstates of the Hamiltonian for the sys- 
tem with M, and thus suitable for orthogonal set of bases 
in the Fock space of polaritons. When the matter config- 
uration transforms to M! , definition of polaritons alters 
accordingly, and the number states of the polaritons in 
M! should be used as new bases. As a result of change 
of bases, the state vector undergoes unitary transforma- 
tion. Thus, when the matter configuration continuously 
transforms as M. (t) and we insist on using the Fock state 
bases in instantaneous matter configuration, the state 
vector undergoes extra evolution in addition to the usual 
time evolution. This representation (the bases follow the 
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eigenstates of the time-dependent Hamiltonian instanta- 
neously) is often employed in the adiabatic approxima- 
tion |38) . The effective Hamiltonian describin g su ch state 
evolution can be derived by a few ways |0, [121 [l4| |29fl 
that are equivalent each other, and here we follow the 
one in Ref. 

The Schrodinger equation is written by 

ih^- t \m) = H(M(t))\m) 

= [H v=0 (M(t)) + K(M{t))]\m), (36) 

where H v= q is the part that does not include the velocity 
of the matter explicitly, and K = H — H v=0 . Expanding 
the state vector as \tp(t)) = J2i c i(^)|^i(^(*))) with the 
eigenvectors satisfying 

[H v=0 (M(t)) - Ei(M(t))]\MM(t))) = 0, (37) 

Eq. ( p6| ) becomes 

ihc^t) = E iCi {t) + y ^{i> i \k\i> j )c j {t) 

3 

i 

= J2(^\H cS \^)cAt), (38) 



where we omit the argument M. (t) for brevity. The sym- 
bol d/dftA indicates differentiation with respect to the 
positions of the matter, e.g., d/dL in the situation of 
Fig. [j] The matrix element of the effective Hamiltonian 
is thus given by 



dependent part 



(^|# eff |^-) = EAj + (il>i\K\il>j) 



-ihM{t){ip t 



d 
dM 



(39) 



Now we express the effective Hamiltonian by using 
the instantaneous creation and annihilation operators 
b^Mit)) and b n (M(t)) (we omit the argument M(t) 
below) . The first term on the right-hand side of Eq. ( |39| ) 
corresponds to hLj r Jr n b n in the effective Hamiltonian. 
Substituting the field expansions (B5T) into the velocity- 



K = J dv{ep(r,t)±(r,t) ■ [v(r,i) x B(r,t)] 



-Y(r,t).[v(r,t)-V]X(r,t)}, 



we obtain 



(40) 



(41) 



where we defined 



F, 



(i) 



2 V LO n 



J dr{ep(r, t) [v(r, t) ■ [X n ,(r) ■ V]A n (r) 



-X„,(r)-[v(r,t)-V]A n (r) 

-Y n (r)-[v(r,f)-V]X n ,(r)}. (42) 

The last term on the right-hand side of Eq. ( ^9|) is 
treated as follows. Differentiating the eigenvalue equa- 
tion (|37|) with respect to M, we find 



E; — Ej 



(43) 



for i 7^ j, and we take the eigenstates as (ipi\ = 0. 
From Eqs. © and ©, a term bj>l, in ^f-, for 

example, corresponds to ib n b n , / (w„/ — u> n ) in the effec- 
tive Hamiltonian, since Ei — Ej = h(uj n ' — uj n ) for non- 
vanishing matrix element. The correspondences of terms 

dH v=0 
dM 

fore obtained as 



(b n ±bi)(b n ,±b~i) 



— Cb^-b 2 ) 



U)„ + LU n i 

i 



±- 



(44a) 
(i>U>l - b n ,b n ) 

-(bib n ,-h n bi). 

(44b) 



In the Hamiltonian (|27]), the density p depends on the 
matter configuration, and then 



dM 



dt 



U(r,t)-ep(r,t)X(r,t) •X(r,t)- 



1 



2mp 2 (r, t) 



Y 2 (r,i) 



mil 2 



X 2 (r,t) . (45) 



Substituting the field expansions ( p5| ) into Eq. ([15]), and applying Eqs. (Q), we find that the last term on the 
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right-hand side of Eq. (|39j) becomes in the effective Hamiltonian as 



+1 



where we defined 



p(2) _ 1„(3) 



^—(6.^1 - U„) - + V 3 j) -J_(St S n - - W-t,) 

+ w n ' \ ^ j u> n — uj n i 



mQ 2 



F« = | dr ^M|_£_ [ n n (r)- e p(r ) t)X n (r)].X n Kr)-^X„(r) -X..,!.-) 



(3) 



r/r 



^(r.i) 



-Y n (r).Y n ,(r). 



^ 2mp 2 (r,t) 

Rearranging the terms in Eqs. ( f4l| ) and (46), we obtain the effective Hamiltonian for polaritons as 



with 



, 1 F (2) _ 1 p (3) 
£ nn i 2 nn 4 nn 

■ r ««' " r .,,2 : 



(n = n') 



,(46) 

(47a) 
(47b) 

(48) 



(49) 



This effective Hamiltonian for the polaritons with mov- 
ing matter is the main result of the present paper, ft 
is interesting to note that the effective Hamiltonian for 
polaritons ( f48| ) and that for photons based on the ex- 
ternal boundary conditions @ have a similar form with 
respect to the creation and annihilation operators, sug- 
gesting that Eq. ( f48|) reduces to Eq. @ in some limiting 
case. This will be explicitly shown in Sec. IV B for one- 



dimensional case. It is also suggested that the squeezed 
state of polaritons will be generated by oscillation of the 
matter at an appropriate frequency by analogy with the 
case of photons j^, [TO] . The important difference be- 
tween Eqs. ( pL8| ) and ([)]) is that the time evolution is fully 
described in a common Hilbert space in Eq. (0) in con- 
trast to Eq. (^|) in which the Hilbert space changes by 
the mirror motion. 



One-dimensional case 



In order to compare our effective Hamiltonian (Mg) 
with Eq. (^]), we consider the one-dimensional moving- 
mirror problem as illustrated in Fig. [l| We assume that 
the system is uniform in the x and y directions, and con- 
sider only the x components of the vector fields A x (z), 
IL x (z), X x (z), and Y x (z) without loss of generality (we 
omit the subscript x below). The normal- mode equations 



( |3l"l ) reduce to 

n„(z,t) 
n„(z,t) 



e A n (z,t) + ep(z,t)X(z,i), (50a) 
' K(z,t), (50b) 



Y n (z,t) = mu 2 n {t)p{z,t)X n {z,t) 1 (50c) 
Y n (z,t) = p(z,t) [mn 2 X n (z,t) + eA n (z,t)] .(50d) 



From these equations, we obtain 



"lit) 



where 



iOM) = 1 



e n (z,t)A n (z,t) = 0, 



2 p(z,t) 1 
e m uJ 2 (t) — O 2 



(51) 



(52) 



can be regarded as the dielectric constant. Thus, in our 
effective Hamiltonian, the dispersion relation is included. 
If the reservoir is taken into account, the Kramers-Kronig 
relations will be satisfied as shown in Ref. |36| for static 
dielectrics. 

The properties of polaritons in the matter significantly 
depend on the sign of the dielectric constant e n . When e n 
is negative, the wave function of the polariton decays in 
the matter, and then polaritons localize between mirrors 
and the energy spectrum is discrete. This condition is 
given by f2 2 < lu 2 < SI 2 + lu 2 , where lu p = (e 2 p/e m) 1 / 2 
is the plasma frequency. When e n is positive, the wave 
function extends indefinitely inside the matter and the 
energy spectrum is continuous. 
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Let us consider the case in which the matter is uni- 
form, i.e., p(z, t) = p6(—z) + p9(z — L(t)), where p is the 
polarization density in the matter and 8(z) is the Heav- 
iside function. In this case, Eq. ( |5l| ) can be solved, and 
when e„ is negative in the matter, the solutions are given 
by 



A n (z,t) = 

a n (t)e Kn ^ z cos 
a n (t) cos k n (t) 



k„(t)L(t) 
2 

_ L(t) 
2 



,-K„(t)[z-L(t)] 



COS 



k n (t)L(i) 
2 



for the even number of nodes, and 
A n (z,t) = 

-a n (f)e K " (t ) z sin 



(f) sin k n (t) 



a n (t)e 



k n (t)L(t) 
2 

2 



-re„(t)[*-Z(t)] 



sm 



k n {t)L{t) 
2 



(z<0) 

(0<z< L(t)) 
(z > L{t)) 

(53a) 



(z<0) 

(0 < z < L{t)) 
(z > L(t)) 

(53b) 



for the odd number of nodes, where k n = to n /c, n n = 
l^n^^kn, and we take the label n to be the number of 
nodes. (The solutions for the continuous spectrum are 
given in Appendix ^].) From the orthonormal relation 
J dz [A n (z)H n > (z) ~ Y n (z)X n i(z)] = —dnn', the normal- 
ization constant a n is obtained by 



£0 



Lit) 



"lit) 



K n [t) U 2 n {t)-U? 



(54) 



where the sign of a n is taken to be A n (0,t) > 0. The 
eigenvalues Lu n (t) are determined so that the solutions 
are smoothly connected at z = and z = L(t), giving 
the eigenvalue equation 



2c 



± 



±1/2 



(55) 



where the signs 



and 



correspond to the solutions 
(pa|) and ( |53b| ), respectively. Using Eqs. @, @, f53|), 
(55j), and dp(z, t)/dt = —pL(t)5(z — L(t)), we obtain the 
coefficient (W9) as 



and 




(56a) 



(56b) 



for n ^ n' , where we omit the argument f for brevity. 

First we consider the case in which the time scale of 
mirror motion is much larger than the inverse of the 
plasma frequency uZ . In this case, the transition be- 
tween the discrete and continuous spectrum can be ig- 
nored, and then the continuous spectrum is irrelevant. In 
order to see the relation between our result and Eq. (||), 
we consider the case of metal, which is obtained by set- 
ting SI — 0. The coefficients (l56j) reduce to 



Cnn(t) 



1 



L(t) 



4L(t) 



(57a) 



At) 



and 



C nn i (f ) 



(-l) n+n '0J n (t)LU n ,(t) W n ,(t) 

<4(t)-a£,(t) W*) 



= (57b) 



for n 7^ n' . When lu u <C lu p , which corresponds to the 
case in which the penetration depth of the EM field is 
much smaller than its wave length, uj n and the coefficients 
( f57| ) can be expanded with respect to 77(f) = c/[L(t)uj p ] <C 
1, giving 



*(*) = 



C nn {t) — 



L(t) 

1 - 
L(t) 



l-2 V (t)+4 V 2 (t) 



24 



v 3 (t) 



0(r/ 4 ), 



(58a) 



[1 - 277(f) +4r7 2 (t) 



C nn i (f ) 



4L(f) 

(8 + nV)r7 3 (f)] +0(rf), 

L(t) (-l) n+n 'nn' fnf 
n 



L(t) n 2 - n' 2 
-i(24 + mV)77 3 (f) 



(58b) 
l-2ry(f)+47j 2 (f) 

-0(?7 4 ). (58c) 



If we identify the photon operators a„ in H° s (t) [Eq. (|g)] 
as the polariton operators b n , we find 

Hf(t) = [1 - 277(f) + 477 2 (f)]i?f (f) + 0(rf). (59) 

When we neglect the terms of order 0(77), the effective 
Hamiltonian for polaritons H£ s (t) reduces to that based 
on the external boundary condition H° (f), and there- 
fore, our method reproduces the existing results of the 
moving-mirror problem in the limit of 77 — > 0. It is in- 
teresting to note that H% (t) is proportional to i?° ff (f) 
up to the second order of 77. This physically indicates 
that the time scale is delayed by the factor ~ 1 — 277 due 
to the coupling of the photon field with the matter field, 
i.e., the EM field drags electrons in the mirrors when 
it is excited. In other words, photons in the cavity are 
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COp/CO c 

FIG. 2: The eigenfrequencies of polaritons ujo/u! c (solid line) 
and u>i /cue (dashed line) and the coefficients of the effective 
Hamiltonian —CooL/L (dotted line) and CoiL/L (dot-dashed 
line) as functions of ui p /u c , where uj c = 2nc/L. 



dressed by plasmons in the cavity mirrors, forming the 
cavity polaritons. 

Figure || shows the lowest two eigenfrequencies deter- 
mined from Eq. ( |55| ) with fl = and the coefficients 
([57]) for n = and n' = 1 as functions of the plasma 
frequency normalized by lu c = 2ttc/L. We find that the 
eigenfrequencies and the coefficients become small when 
ujp is comparable to cu ci and they go to the asymptotic 
values for lu p 3> uj c . For example, in the case of a su- 
perconducting microwave cavity, the plasma frequency is 
in the X-ray region and io c is in the microwave region, 
and therefore the correction to the ideal metal can be 
neglected. However, this correction might be significant 
if Lu p becomes small due to, e.g., decrease of the carrier 
density. 

When the time scale of mirror motion is comparable 
to oj ~ , transition between discrete and continuous spec- 
trum occurs. The transition is significant when the mir- 
ror vibrates at the frequency ujm > uj p — ui n , which results 
in decay of the polaritons in the nth mode into contin- 
uous spectrum, namely, photons leak out of the cavity. 
When the mirror moves as L{t) = Lq + isimuMt with 
f < Lo, the time-dependent part of the effective Hamil- 
tonian can be written by V e liA)Mt + V^e~ lUMt , where V 
is obtained by replacing L(i) with £u>m/2 in H^ s . Using 
Fermi's golden rule, the decay rate of a photon in the nth 
mode is estimated to be 



show that when Q = and w n < w p < the decay 
rate ( |6C| ) reduces to R n {uj M ) ~ n£ 2 ujp/(4LQU! M ). 

In the effective Hamiltonian H° s (t), all the energy lev- 
els are commensurate (uj n = mr/Lo), and the created 
photons make transition to higher levels unlimitedly as 
Lu n — > U2n — > • • • due to the resonance. In our one- 
dimensional effective Hamiltonian H^ s (t) [Eq. (|59|)], on 



the other hand, the transition stops at 



due to the 



incommensurate energy levels, or the decay into the con- 
tinuous spectrum occurs, and thus the resonant enhance 
ment of the DCE is to be naturally suppressed at ~ lo 



2 

Ctl i=l,2 



W|2 



kn 



(60) 



V. CONCLUSIONS 

We formulated the DCE in terms of microscopic field- 
mattcr theory, in which the EM field and the polarization 
field in the matter are treated on an equal footing. This 
enabled us to study the DCE without boundary condi- 
tions and without changing the Hilbert space. We de- 
rived the effective Hamiltonian for polaritons with mov- 
ing matter, and applied it to the one-dimensional cavity 
with a moving mirror. We obtained the corrections to the 
results based on the external boundary conditions: the 
time scale of the dynamics is delayed when the plasma 
frequency is comparable to the resonant frequency of the 
cavity. This effect is attributed to the fact that the pho- 
tons in the cavity are dressed by the electrons in the 
mirrors. 

Finally, we comment on possibility of experimental ob- 
servation of the DCE. This effect has not been demon- 
strated in laboratories yet, since time scale of the phe- 
nomena is extremely fast. An efficient way to observe 
the DCE is to accumulate photons in a cavity by vibrat- 
ing the mirror surface of the cavity at twice the reso- 
nant frequency PQ |. However, the resonant frequency of 
the high-Q cavity is typically > 10 GHz, and it is quite 
difficult to excite oscillation at such a high frequency. 
One possibility to overcome this obstacle might be to 
slow down the speed of light. Using the electromagnet- 
ically induced transparency in an ultracold atomic gas, 
the speed of light can be reduced to ~ 10 m/s in the 
regime of visible light jf(| . If this technique can be ap- 
plied to much lower frequency, the resonant frequency of 
the cavity is significantly reduced, which enables us to 
resonantly vibrate the cavity wall to observe the DCE. 



APPENDIX A: DERIVATION OF THE 
EFFECTIVE HAMILTONIAN (S) 



where U fc is the matrix element of V with respect to 



the nth mode in the discrete spectrum ( J53| ) and the mode 
labeled by k = (uj n +ujM)/c and i = 1 , 2 in the continuous 



spectrum (CI). Here we took the wave number k between 
mirrors as the mode index in the continuous spectrum. 
The explicit form of vj^ is given in Appendix |C[ We can 



In order to make this paper self-contained, we de- 
rive the effective Hamiltonian (^|) for the DCE based on 
the external boundary condition following the method 
in Ref. p^ |. We restrict ourselves to the moving-mirror 
problem, while in Ref. [jl2| change of the dielectric con- 
stant is also considered. 
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The equations of motion are given by and the boundary conditions are A(0, t) = A(L(t),t) = 0. 

The field operators are expanded as 

-A(z,t) = -E(z,t), (Ala) 



dt 

8 - 3 2 



^-E(z,t) = -c 2 —A(z,t), (Alb) 



A(z,t) = J2 

n 

E(z,t) = J2 



2e u! n (t) 



(f> n (z,t)[a n (t) + al(t)}, 



2e 



<f> n (z,t)[a n (t) ~ al(t)}, 



(A2a) 
(A2b) 



where uj n (t) = nitc/ L{t) and 



4>n{z,t) 



■ sin 



Lit) \L{t) 



(A3) 



From Eqs. (A2) and ([A3]), the annihilation operator is written by 
a n {t) - 



L(t) 1 

dz-<f> n (z,t) 



2e uj n (t) 



A{z,t)-u 



2e 



hu} n (t) 



E(z,t) 



Using Eqs. (Al), time derivative of Eq. (A4) reads 

da n (t) 
dt 

L{t) 1 







dz-(f> n (z,t) 



2 -^^E(z > t) + u 



2ea 2 d 



huj n {t) dz 



c^A(z,t) 



(A4) 



L ^ f m A X A ( A / 2 W*) % , 



L(t)nir f L{t) , 1 
dz—. 



L 2 {t) J 2 V L{t) ' 



L(t) 



2 J^A(z,t)- l 



2e 



•E{z,t) 



L(t), u . L(t) ^ (-!)»+" W 



TlUJ n {t) 



(A5) 



Thus we find that the effective Hamiltonian (^|) gives this On the other hand, substituting Eqs. (|jl|a) and (|3l|d) 
time evolution ( |A5| ) by the Heisenberg equation. into A n / and Y n > in , and using V • II„ = 0, we find 

I 11 fi ' — Ifh ' 7h ' \ ^ ) 

APPENDIX B: ORTHOGONALITY OF THE 
MODE FUNCTIONS 

We give a proof of the orthogonality of the mode func- 
tions in Eq. (|32|). We consider the integral 

I nn , = J dr[A n ,(r)-U n (r)-Y n ,(r)-X n (r)}. (Bl) 



i.C, Inn' 



From Eqs. (B2) and (B3h, we obtain {w* - u 2 n ,)I n 



for u>* ^ uj^, 



APPENDIX C: CALCULATIONS FOR 
CONTINUOUS SPECTRUM IN 
ONE-DIMENSION 



Substituting Eqs. (pip) and (Bib) into II ra and X„ in 

2 



u*I nn >, we find 



When s n is positive in the matter, the energy spectrum 
is continuous, and we use the wave number k between 

(B2) 

mirrors as the mode index. There are two independent 
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solutions of Eq. ( f5l| ) for given k as 



4H*) 



4 2) (*) 



ockyz cos ^ sin nz — sin ^ cos kz) (z < 0) 

o; fc sinfc(z - L/2) (0 < z < L) 

a k [| cos sin k(z - i) + sin ^ cos k(z - L)] (z > L), 

f3 k (~ sin ^ sin kz + cos ^ cos kz) (z < 0) 

(0 < z < L) 
L) + cos ^j- cosk(z - L)] (z>L), 



/3 k cosk(z- L/2) 
(3 k [-|sin^sinK(z- 



(Cla) 
(Clb) 



1/2 

where k — ej k, and the normalization constants 



a k 



(3 k = 



1 9 kL 9 kL 
7reoV £ fe ( — cos — + sin — 

e k 2 2 



1 9 fc-Zy 9 kL 

7T£oV £ fc | — sin — + cos — 
e k 2 2 



-1/2 



(C2a) 



-1/2 



(C2b) 



are determined from the orthonormal relation 

dz[A®W$(z)-Y®(z)xip{z)] 



8{k — k')6ij . 

(C3) 



We can show that the functions (p3[) and (CI) are orthog- 



onal each other. The matrix element in Eq. ( |60| ) is 
obtained by the straightforward calculation as 



V, 



,(i) _ i ^ n ihe £ojM 

kn ~ \ 1 ) a 



k 

4 V k n 
a k \a n \k p (k + k n ) 



( K 2 + 4)(fc2_ fc 2 )(fc 2_ A .2 ) l/2 
9 /cZ/ 9 . \ 

k K n cos — kk^ sm — I 



2 ; ' 



(C4) 



where fco = Q/ c , k p = u) p /c, and V^ n is obtained by 
the replacement a k — > j3 kl cosfcL/2 — > — sin/cL/2, and 
sin/cL/2 — * coskL/2. 



[1] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 

(1948). For review see, for example, G. Plunien, B. 

Miiller, and W. Greiner, Phys. Rep. 134, 87 (1986). 
[2] L. Parker, Phys. Rev. Lett. 21, 562 (1968); Phys. Rev. 

183, 1057 (1969). 
[3] G. T. Moore, J. Math. Phys. 11, 2679 (1970). 
[4] S. A. Fulling and P. C. W. Davies, Proc. R. Soc. London 

A 348, 393 (1976). 
[5] See, for example, M. Castagnino and R. Ferraro, Ann. 

Phys. (N.Y.) 154, 1 (1984). 
[6] V. V. Dodonov, A. B. Klimov, and V. I. Man'ko, Phys. 

Lett. A 142, 511 (1989). 
[7] M. T. Jaekel and S. Reynaud, J. Phys. I France 2, 149 

(1992) . 

[8] P. A. Maia Neto and S. Reynaud, Phys. Rev. A 47, 1639 

(1993) ; P. A. Maia Neto, J. Phys. A 27, 2167 (1994). 

[9] V. V. Dodonov, A. B. Klimov, and V. I. Man'ko, Phys. 

Lett. A 149, 225 (1990). 
[10] S. Sarkar, Quantum Opt. 4, 277 (1992). 
[11] M. Razavy and J. Terning, Phys. Rev. D 31, 307 (1985). 
[12] C. K. Law, Phys. Rev. A 49, 433 (1994). 
[13] C. K. Law, Phys. Rev. A 51, 2537 (1995). 
[14] R. Schiitzhold, G. Plunien, and G. Soff, Phys. Rev. A 57, 

2311 (1998). 

[15] E. Yablonovitch, Phys. Rev. Lett. 62, 1742 (1989). 
[16] V. V. Dodonov, A. B. Klimov, and D. E. Nikonov, Phys. 
Rev. A 47, 4422 (1993). 



[17] T. Okushima and A. Shimizu, Jpn. J. Appl. Phys. 34, 
4508 (1995); A. Shimizu, T. Okushima, and K. Koshino, 
Materials Science and Engineering B48, 66 (1997). 

[18] H. Saito and H. Hyuga, J. Phys. Soc. Jpn. 65, 1139 
(1996); ibid. 65 3513 (1996). 

[19] G. Barton and C. Eberlein, Ann. Phys. (N.Y.) 227, 222 
(1993). 

[20] G. M. Salamone and G. Barton, Phys. Rev. A 51, 3506 

(1995) . 

[21] G. Barton and A. Calogeracos, Ann. Phys. (N.Y.) 238, 
227 (1995); 

[22] G. Barton and C. A. North, Ann. Phys. (N.Y.) 252, 72 

(1996) . 

[23] R. Giitig and C. Eberlein, J. Phys. A 31, 6819 (1998). 

[24] C. Eberlein, J. Phys. A 32, 2583 (1999). 

[25] V. V. Dodonov, A. B. Klimov, and D. E. Nikonov, J. 

Math. Phys. 34, 2742 (1993). 
[26] A. Lambrecht, M. T. Jaekel, and S. Reynaud, Europhys. 

Lett. 43, 147 (1998). 
[27] G. Plunien, R. Schiitzhold, and G. Soff, Phys. Rev. Lett. 

84, 1882 (2000). 
[28] J. Schwinger, Proc. Natl. Acad. Sci. USA 89, 4091 

(1992); 90, 958 (1993); 90, 2105 (1993); 90, 4505 (1993); 

90, 7285 (1993); 91, 6473 (1994). 
[29] C. Eberlein, Phys. Rev. Lett. 76, 3842 (1996); Phys. Rev. 

A 53, 2772 (1996). 
[30] V. V. Dodonov and A. B. Klimov, Phys. Rev. A 53, 2664 



12 



[31] 
[32] 
[33] 
[34] 
[35] 

[36] 



quant-ph/0106081 



(1996). 

For recent review, V. V. Dodonov 

M. Koashi and M. Ueda, Phys. Rev. A58, 2699 (1998). 
U. Fano, Phys. Rev. 103, 1202 (1956). 
J. J. Hopfield, Phys. Rev. 112, 1555 (1958). 
R. J. Glauber and M. Lewenstein, Phys. Rev. A 43, 467 
(1991). 

B. Huttner, J. J. Baumberg, and S. M. Barnett, Euro- 
phys. Lett. 16, 177 (1991); B. Huttner and S. M. Barnett, 
Europhys. Lett. 18, 487 (1992); Phys. Rev. A 46, 4306 



[37] 

[38] 
[39] 



[40] 



(1992). 

R. Matloob, R. Loudon, S. M. Barnett, and J. Jeffers, 
Phys. Rev. A 52, 4823 (1995). 

A. Messiah, Mecamque Quantique (Dunod, Paris, 1964). 
L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Elec- 
trodynamics of Continuous Media (Pergamon, Oxford, 
1984). 

L. V. Hau, S. E. Harris, Z. Dutton, and C. H. Behroozi, 
Nature 397, 594 (1999). 



